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Abstract 

In this review paper the QCD vacuum properties and the struc- 
ture of color fields in hadrons are studied using complete set of gauge- 
invariant correlators of the gluon fields. Confinement in QCD is pro- 
duced by the correlators of some certain Lorentz structure, which vi- 
olate abelian Bianchi identities and therefore are absent in the case of 
QED. These correlators are used to define an effective colorless field, 
which satisfies Maxwell equation with nonzero effective magnetic cur- 
rent. With the help of the effective field and correlators it is demon- 
strated that quarks are confined due to effective magnetic currents, 
squeezing gluonic fields into a string, in agreement with the "dual 
Meissner effect" . Distribution of effective gluonic fields are plotted in 
mesons, baryons and glueballs with static sources. 



1 Introduction 



The QCD is a unique example of field theory, lacking internal contradictions 
and at the same time explaining all physical phenomena in strong interactions 
[H 12]. The theoretical understanding of QCD is complicated due to the 
fact that all its basic features are of the nonperturbative nature, and the 
QCD vacuum is a dense and highly nontrivial substance. In fact, in the 
modern quantum field theory one often represents the vacuum as a specific 
material substance with definite characteristics in direct analogy with the 
condensed matter physics. As illustrating examples one can mention the 
Casimir effect and relative phenomena, and also the Higgs mechanism in 
the standard model. In the last case one deals with the vacuum condensate 
of the scalar field (0), while quantum excitations above this condensate are 
considered as Higgs particles. 

The nontriviality of the QCD vacuum is revealed by the fact that this 
medium has nonzero values of gluonic condensate [Sj, (Fl^u-^li^) = (600MeV)'', 
and of the quark condensate, (qq) = —(250 MeV)^. As it has become clear 
during last decades it is the vacuum properties which bring about confine- 
ment (see, e.g., the review For theoretical calculations in QCD one usu- 
ally exploited till recently the perturbation theory augmented by some models 
of nonperturbative mechanisms. The situation changed with the advent of 
the QCD sum rule method [3j, which uses the gauge-invariant formalism of 
condensates to describe the nonperturbative contributions. However for most 
effects at large distances this method is not sufficient, e.g. for confinement or 
spontaneous violation of chiral symmetry. The systematic description of all, 
in principle, QCD phenomena is made possible due to the appearence of the 
Vacuum Correlator Method (VCM), see IHl El and the review 0, which 
exploits as basic elements the complete set of field correlators of the form 

^Si.../.„;.„(a;i,...,a;„,a;o) = {Ti G^,^,{xi,xo)...G^„^„{xn,Xo)) (1) 

where the notation G^^^^j^^Xi, Xq) is used for the gluonic field strength covari- 
antly shifted along some curve, see Eq. (0). 

The basis of the VCM are the gauge-invariant Green's functions of white 
objects, which can be written as path integrals through field correlators (jT]) 
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using the cluster expansion, see e.g. [7| El CDl- A question may arise at 
this point, why one considers in VCM only white (i.e. gauge-invariant) ob- 
jects, and not for example propagators in some fixed gauge? The answer is 
tightly connected to the difference between gauge invariance in abelian and 
nonabelian theories. In the abelian theory, e.g. in QED, the requirement 
of gauge-invariance does not forbid to consider the problems with formally 
gauge-noninvariant asymptotic states, like electron-electron scattering. The 
gauge invariance of the cross-section occurs in this case due to the conser- 
vation of the abelian current. In the nonabelian theory with confinement, 
like QCD, the situation is different and the problem of scattering of isolated 
quarks has no sense. Formally one can see that nonlocal gauge noninvariant 
matrix element vanishes when being averaged over gluonic vacuum, which is 
a property of the nonabelian SU{N) group. Therefore instead one consid- 
ers in QCD the quark-quark scattering for quarks inside white objects (i.e. 
described by the gauge- invariant functions), such as hadrons. The same is 
true for the problems connected to the spectrum of bound states — while 
in QED the problem of a neutral atom spectrum is just as valid as that of 
the spectrum of a charged ion, in QCD an analog of the last problem has no 
meaning. 

Therefore the set of correlators (0) can be considered as a starting dynam- 
ical basis yielding a phenomenological gauge-invariant description of physical 
processes. In fact, however, the situation is much more interesting. First of 
all, the lattice calculations give important evidence that already the first 
nontrivial correlator with n = 2 dominates, and the total contribution of all 
highest correlators is below few percent, see review |11J. As was shown in 
jni Ej ; the lowest (or Gaussian^ as it will be called in what follows) correlator 
can be expressed through two scalar formfactors D{xi — X2) and -Di(xi — X2). 
Secondly, both formfactors have been measured in the lattice calculations 
and have nonperturbative parts of exponential shape with a characteristic 
small correlation length A. Finally, the function D^'^^ (and therefore also 
D{xi — X2), Di{xi — X2)) are directly connected to the Green's functions of 
the so-called gluelumps [121 CHI HI] • The latter can be calculated analytically 

^using the analogy with so-called Gaussian, or white, noise described by a quadratic 
correlator, in this case with vanishing correlation length 



3 



in terms of the only mass scale of QCD, e.g. through the string tension a 
and the coupling constant a^. Thus the formulation of the nonperturbative 
dynamics in QCD turns out to be selfconsistent and one should in addition 
calculate A through a, what was done earlier in ^3], and also connect a and 
^QCD and write down the explicit form of correlators D^"'\ This way one 
would also be able to understand analytically the dominance of D^'^^ (one 
can find the first results in this direction in [TBI I17j ). 

With all that the formalism of field correlators is to a large extent unusual 
to physicists, brought up in the standard lore of perturbative, or even more, of 
abelian gauge theory. In the context of the confinement problem such a "lin- 
ear" abelian approach is realized in the so-called "dual Meissner scenario", 
which contains a simple qualitative picture of the confinement mechanism 
in QCD jlSl ITU] . In this approach the acting roles have charges (quarks) 
and the monopole medium filling the vacuum. Many lattice and analytic 
studies, see, e.g., [20] - |22]; demonstrate that the string formation between 
quark and antiquark is connected in this picture with the appearence of cir- 
cular monopole currents k around the string, which obey the dual Ampere 
law k = rot E. From the physical point of view this situation is similar to 
the Meissner effect in the standard superconductivity phenomenon, modulo 
interchange of effective electric and magnetic charges. On the other hand, 
the defect of this picture is that the very notion of the magnetic monopole 
cannot be exactly defined in QCD. This arbitrariness can be seen, first of 
all, in the gauge dependence of the monopole definition, and secondly, in the 
difficulties with the continuum limit for the lattice monopoles, defined by the 
flux through an elementary cube. There is a lot of literature, with different 
suggestions on how to deal with this problems, see e.g. |24j . 

While confinement properties are studied on the lattice numerically, in- 
cluding partly the abehan projection, they are also an object of investigation 
in the effective Lagrangian approach and in different dielectric vacuum mod- 
els of QCD |2n]-jSll- The basic field theory problem in this case is replaced 
by a classical variational problem for the effective Lagrangian, which yields 
a system of differential equations, to be solved numerically. In this way one 
introduces an effective dielectric constant of the vacuum, depending on the 
effective fields and ensuring quark confinement. 
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In what follows we shall use another approach which is fully gauge- 
invariant and yields a simple and selfconsistent picture of the confining string 
formation. Namely, using the field correlator method as a universal language, 
one can define gauge-invariant (with respect to the gauge symmetry of the 
original nonabelian theory) effective field J-'/iuix) via the W-loop. The 
effective electric field near the charge turns out to be the gradient of the 
color-Coulomb field, and in the case of an abelian theory J-'^u{x) is the stan- 
dard field strength. The effective field satisfies Maxwell equations, having on 
the r.h.s. electric current and magnetic current k^. The source of A;^ is 
primarily the triple correlator of the form (EEB) (as was already found in 
0]) describing the emission of the color-magnetic field by the color-electric; 
the latter can be visualized as the emission of the color-magnetic field by an 
effective magnetic charge (monopole). In the language of field correlators one 
can easily demonstrate that the system of equations for the effective fields 
describes the QCD string and the circular magnetic currents around it. In 
this way the picture of the dual Meissner effect is given in gauge-invariant 
terms. 

With the help of one can investigate in detail the structure of the 
QCD string. The first computations of the string profile in have demon- 
strated a very good agreement of the results calculated via D(x^),Di(x^) 
and those obtained independently on the lattice. The following study of the 
string structure j2Hl has shown an interesting phenomenon of the profile sat- 
uration, where the profile (i.e. the field distribution across the string) does 
not change for long enough strings. The relief of the baryon field has turned 
out to be even more interesting. Baryons, and more exact, nucleons are the 
basis of the bulk of the stable matter around us. The physical problem of 
the structure of the baryon field is especially interesting both from theoret- 
ical and practical points of view. Two types of baryon field configurations 
were discussed in the literature: with the string junction in the middle (the 
F-shape) and of the triangular shape (the (A-shape). Using vacuum corre- 
lator method the baryon configuration was computed analytically in [211 110] , 
where the presence of the string junction in the field distribution was ex- 
plicitly demonstrated, thereby excluding the A-type configuration. On the 
other hand, the latter is possible for the three-gluon glueballs and the corre- 
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spending field was calculated in jlU]. One should mention that these baryon 
field distributions are also in agreement with the lattice calculations using 
the abelian projected QCD jH], see also the review paper of Bornyakov et. 
al. 021. 

The field sources in three-gluon glueballs are three valence gluons. The 
field structure of these systems has some specific features, and can be of both 
types, of the A- type (unlike baryons), and of the F-type (like baryons), and 
its study helps to understand better physics of confinement. Moreover, the 
three-gluon glueballs have to do with the processes of the odderon exchange 
(i.e. glueball exchange with odd charge parity), and hence are also interesting 
from the experimental point of view. Therefore, in addition to the effective 
field distributions, we shall also discuss below the W-loops and the static 
potentials of baryons and three-gluon glueballs. 

The paper has the following structure. In chapter 2 the discussion of 
field correlator properties in QCD is given, and in particular the important 
phenomenon of the Casimir scaling is explained. In chapter 3 the effective 
field and currents j^, are introduced and the dual Meissner effect is 
demonstrated. In chapter 4 the static potentials and field distributions in 
baryons and three-gluon glueballs are given. In Conclusions the main results 
are summarized and some prospectives are outlined. 

Everywhere in what follows, if it is not especially stressed otherwise, the 
Euclidean metrics is used with notations for 4-vectors k = (fci, k2, ks, ^4) and 
for scalar products kp = kfj^p'^d^. The three-dimensional vectors are denoted 
as k = {ki,k2,ks) and the Wick rotation corresponds to the replacement 
/c4 — iko- 

2 Properties of QCD vacuum in gauge- invariant 
approach 

2.1 Definition of gauge- invariant correlators 

The following remark is to be made before we proceed. There is an important 
difference between pure Yang-Mills theory (gluodynamics) and QCD, namely 
the latter contains dynamical fermions, in particular light u and d quarks. 
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This circumstance plays no crucial role in the description of confinement since 
gluodynamics confines color as QCD does, which is supported by direct lattice 
calculations (see, e.g. PH]) and different qualitative arguments. Because of 
that in most cases we consider pure Yang-Mills theory in this review, while 
quarks play a role of external sources. 

One of the main objects in gauge theory is the Wegner- Wilson loop ^Sl HE] 
which we denote here as W-loop: 

W{C) = P exp ig j Al{z)t''dz^, (2) 
o 

where t"" - generators in the given representations of the gauge group. W- 
loop defines external current J which corresponds to a point particle charged 
according to the chosen representation and moving along the closed contour 
C. Phase factor for non-closed curve connecting points x and y is also of 
importance 

y 

y) = V exp ig j Al{z)edz^ (3) 

X 

Under the gauge transformations we have 

y) ^ y) = U\x)^x- y)U{y) (4) 

It means that the trace Tr W{C) is gauge-invariant.^ We normalize Tr 
everywhere as Tr 1^ = 1 for the given representation of dimension d. Making 
use of the definition Q, let us introduce G^y^x.xo) as 

Gf,u{x, xo) = $(xo; Xq) (5) 

where F^^ = d^A^ — dyA^ — i[A^Ay\ is nonabelian field strength and the 
curve connecting the points x and xq does not self-intersect. In abelian 
theory G^y{x,xo) = F^y{x), however in Yang-Mills theory G^,^(x, xq) and 
Ffj,u{x) transform differently under gauge transformations, as it is clear from 
(jlj). We can now construct vacuum averages of the products of Gf^^{xn,Xo) 
in the following way 

^ilpa(^) xo) = (Tr G^^{x, xo)Gp^{y, xq)) (6) 



^In the literature the trace is often included in the definition of the W-loop. 
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Djpaapi.^^ 2;, Xo) = (Tr G^u{x, xo)Gpa{y, xo)Gap{z, Xo)) (7) 

and analogously for higher orders. The correlators (jHI), are gauge- invariant, 
but nonlocal - expressions (jH)), © depend on the position of the points x, y, 
z as well as on the position of the point Xq and contour profile used in (0). 
Physical observables such as static potential extracted from the W-loop do 
not depend on xq and contour profiles when all correlators n > 2 are 

taken into account. It is not true, however if one takes only the lowest n = 2 
term. In this case it is convenient to minimize the corresponding dependence 
like one does in perturbation theory minimizing the contribution of omitted 
terms by the proper choice of subtraction point fi on which the exact answer 
should not depend. 



2.2 Computation of the W-loop and Green's functions 
in terms of correlators 

Speaking in general terms, for a given gauge theory each function D^""^ is 
important characteristics of its vacuum structure by itself. What is more 
important, however is the possibility to express W-loop average in terms 
of correlators ©, ©• Indeed, Stokes theorem (or, more precisely, its non- 
abelian generalization |17j-|S2]) leads to 

{TtW{C)) = (riVexpig J xo)f/a^.(^)^ = exp f^(z(7)"A(")[^] 

(8) 

there we have used cluser expansion to exponentiate the series (see, e.g. 
iniEll)- Integral moments A [5] over the surface S of irreducible corre- 
lators, known as cumulants in statistical physics can be expressed as linear 
combinations of the integrals of correlators D^'^\ For example, we have for 
two-point correlator 

A(2)[S] = da,M) J dap„{z2)D^X<^{z,,Z2,x,) (9) 
s s 

For higher terms the ordering is important, see, e.g. flBj, where exact com- 
putations for n = 4 are performed. 
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Expression (jSI) is of central importance for the discussed formalism. Let 
us consider the propagation of the spinless particle with mass m, carrying 
fundamental color charge ("quark") in the field of infinitely heavy "anti- 
quark" jTliniCni- The corresponding gauge-invariant Green's function reads 
as 

g{x,y) = {<p\x)^{x-y)cP{y)) (10) 

where we denote quark field as 0(a;). One can demonstrate that Q{x,y) has 
the following Feynman-Schwinger representation 

g{x,y) = jds J Vz,expi-mh-^jdT(^^^^ • (Tr W^(C)) 

z^,{0)=x^, \ / 

(11) 

where the closed contour C is formed by the quark trajectory z^{t) and that 
of antiquark (the latter is nothing but the straight line connecting the points x 
and y). We have taken spinless case here as the simplest illustrative example, 
for real physical problems with spinor quark fields there is a systematic way 
of analysis of spin effects [SHI IH CHI- The problem of two-body meson state 
or three-body baryon one can be addressed in completely analogous way. In 
all cases Green's function containing full information about mass spectrum 
and wave functions of the system can be re-written in terms of path integrals 
of the W- loops there the latter are expressed via correlators as in (jSI). 

Therefore the set of correlators Z)*^"^ provides rich and, what is more 
important, universal dynamical information one can use to compute different 
nonperturbative effects.^ Let us stress once again that the correlator © 
is itself related to the Green's function of gluon excitation in the field of 
infinitely heavy adjoint source - known in the literature as gluelump 



Coming to practical side of the problem, it is natural to ask what the 
actual behavior of the correlators (jHI), O is and how information about it 
can be gained. This question is simple to answer in perturbation theory 

■^Discussed formalism can be applied in perturbation theory as well. In this context 
it allows to sum up perturbative subseries with subsequent exponentiation with the well- 
known "Sudakov formfactor" as a result of the first approximation, see [HO], 50] and 
references therein. 
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since each D^"-^ is given by perturbative series, see, e.g. [SSlinZl- There are 
a few ways to proceed beyond perturbation theory. The first one is to find 
nonperturbative solutions to the so called BBGKI equations, relating the cor- 
relators of different orders [HI] . This way has brought no essential progress up 
to now. Another analytic strategy suggests to compute correlators in terms 
of gluelump Green's functions J21-Cni- The third and the most successful 
way is to study the problem on the lattice. There are quite a few sets of 
numerical data 1^2]- [HZI, which we discuss below. However it is obvious that 
numerical results concerning one or a few particular correlators are useless if 
general properties of the whole ensemble are unknown. To discuss them we 
come back to the expression (jS)). 

2.3 Gaussian dominance 

It has already been stressed that the price we have payed for manifest 
gauge-invariance of (jH)) is the dependence of (0), (0) on the contour pro- 
files entering These contours are, generally speaking, arbitrary non- 
selfintersecting curves or, better to say, they can be freely chosen in some 
(large enough) set. As a result the quantities A*^")[S'] in (jH)) depend on this 
choice while W{C) is obviously independent on S. The contradiction is spu- 
rious and one can demonstrate that this contour dependence is cancelled in 
the total sum, despite it is present in each individual summand A*^"''[S']. In 
this sense the choice of the surface S* in (jHI) (corresponding to the choice of 
integration contours in the correlators D^'^^) is free, as it should be. We can 
take a different attitude and ask the following question: what is the hierarchy 
of cumulants A^"^[S'] on some particular surface? This question is of general 
interest but it has also important practical meaning - in many problems one 
has to deal with the surface, which is singled out by some physical reasons. 
For a single W-loop it is obviously given by minimal area surface, bounded 
by the contour. In more complicated case of interacting loops [HH] the surface 
corresponding to the minimal energy of the system can be taken. In any case, 
it is instructive to make a distinction between two different scenarios: 

oo 

A(2)[5] > ^ A(")[^] (12) 

n=3 
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which is referred to as stochastic scenario, while the case when (I12|) does not 
hold (for example, all cumulants are of the same order) is known as coher- 
ent. General framework described in the present paper takes into account 
effects of all cumulants but as it should be clear, it shows its strong sides in 
stochastic case. The lowest two-point Gaussian cumulant (jHl) is dominant in 
stochastic ensemble, while higher order terms can be considered as small cor- 
rections. This situation is known as Gaussian dominance. Then we can ask is 
the QCD vacuum stochastic or coherent? To answer this question in straight- 
forward way one has to compute (for example, numerically on the lattice) 
different cumulants and check them against (|12j). Unfortunately this research 
program is too intricate for modern lattice technologies and almost all ac- 
tual results are obtained for Gaussian cumulant only. There are important 
indirect evidences however supporting the idea that Yang -Mills vacuum is 
indeed stochastic and not coherent in the sense of ()12|). Of prime importance 
in this context is Casimir scaling phenomenon [021- [ZS; see also [7S]-[7n]. Us- 
ing (jSI) and taking into account well known relation between static potential 
and W-loop average, one can get, assuming Gaussian dominance 

V{R) = lim ^ g^A^^^ [S = R x T] (13) 

and, according to (jHl) and we have V{R) ~ Cd, where an eigenvalue 
of Casimir operator in the representation d is given by dabf^t^ = Cd ■ Id- 
Let us remind that representation of the Lie group SU{N) of dimension d is 
characterized by A^^ — 1 generators t", which can be realized a.s dxd matrices 
commuting as [tH''] = if'^^'H^. Proportionality of the static potential to Cd is 
called Casimir scaling [77] and was discussed for the first time in [75] . 

It can easily be shown that contributions from higher cumulants to the 
static potential ([T^ are, generally speaking, not proportional to Cd (despite 
they can contain linear in Cd terms). Therefore a good accuracy (devia- 
tion not exceeding 5%) of Casimir scaling demonstrated on the lattice is a 
serious argument in favor of Gaussian dominance. Moreover, attempts to 
reproduce Casimir scaling in many other models of nonperturbative QCD 
vacuum encounter difficulties [73 jEDl HI] ■ Another argument is the observed 
independence of radius of the confining string between quarks on their non- 
abelian charge (i.e. on representation d) [HT]. These results would look as 
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fine tuning effects without Gaussian dominance. It is also worth mentioning 
that "vacuum state dominance" successfully used for years in QCD sum rules 
formalism is nothing but Gaussian dominance in our language. 

2.4 Structure of two-point correlators 

We have mentioned above the relation between correlators D^""^ and gluelump 
Green's functions. For the simplest Gaussian correlator (jH)) this can be seen 
clearly if the contours are straight lines and points x, y, Xq belong to one and 
the same line. The correlator depends on the only variable z = x — y in this 
case and can be represented as 




Expression (jl4p contains phase factor in the adjoint (compare with the pre- 
vious formulas where we worked with fundamental phase factors, i.e. with 
N X N matrices) which makes its physical content self-evident. Namely, 
gauge-invariant function D^'^\z) describes gluon propagation in the field of 
infinitely heavy adjoint charge at the origin in full analogy with fundamental 
case (compare (fTUI) and (HH)). 

Confining string worldsheet given by the surface 5* in (jH)) interacts with 
itself by gluelump exchanges. This interaction depends on the profile of 5* in 
such a way that the total answer for the W-loop average is S*- independent. 
Gaussian dominance means, qualitatively, that for some particular surface 
this "gluelump gas" becomes "ideal" and integral contribution of higher cu- 
mulants A*^"\ n > 2 is small on this surface. This also means that two-gluon 
gluelumps weakly interact with each other. The deviation from Casimir scal- 
ing (as we have already noticed, it is small) can be expressed in terms of 
irreducible averages of gauge- invariant operators (Tr OiTr O2), describing 
interaction of gluelumps ^H]- One immediately realises that such deviation 
is suppressed in large limit. To avoid misunderstanding let us stress that 
gluelumps do not exist as physical particles in the spectrum of the theory. 
It would also be wrong to interpret (fT^ in terms of " massive gluon" . In a 
limited sense gluelumps are analogous to Kalb-Ramond fields which describe 
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dual vector bosons and play important role in constructing string represen- 
tation of compact QED jHS] and abelian Higgs model [82 (see also (S^ ISij). 
The discussed picture with gluelump ensemble on the worldsheet makes sense 
only in the presence of external current, forming the W-loop. On the other 
hand the correlator (|Tl|l may be studied as it is, with no reference to any ex- 
ternal source. Before we discuss actual lattice results, it is useful to represent 
(|14p in terms of two invariant formfactors D{z^) and Di^z^) 0-0 

I f d d \ 

Confinement (linear potential between static quark and antiquark in fun- 
damental representation) takes place in Gaussian dominance picture when 
D{z'^) is nonzero. At large distances we have from (jH)), p3|) for static poten- 
tial V{R) and string tension a: 

V{R) = aR + 0{R^) ■ (^ = \j d^zD{z^) (16) 

while at small distances perturbative contribution dominates j56| 15 7j. Non- 
perturbative part of the correlator is usually taken as 

D(22) ~ exp(-|^|/A) (17) 

and this exponential fit is in very good agreement with lattice data at large 
enough distances. The situation with nonperturbative component of the 
function Di{z^) is less clear. In any case, physically the exact function Di{z'^) 
containing perturbative and nonperturbative pieces must be exponentially 
suppressed at large enough distances. It is important that from practical 
point of view one has no need to know the detailed profile of formfactors 
D[z^), Di{z'^): physical quantities are given as integral moments of these 
functions as in (jl6|) . Quantity A is known as correlation length of QCD 
vacuum and as it is clear from our discussion this quantity is nothing but the 
inverse mass of the lowest gluelump: A = 1/M. On the other hand, typical 
size of vacuum domain where fields are correlated is given by the same A jHH] • 
We use numerical value A = 0.2 fm in accordance with the lattice results. 
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Physics of nonlocality switches on at distances larger that A and has many 
phenomenological manifestations. One of the most interesting - the confining 
string formation - will be discussed in what follows. 

So far we have not mentioned the problem of deconfinement. There are 
basically two groups of physically interesting questions related to this prob- 
lem. The first one covers dynamical aspects of the phase transition, while 
the second group deals with symmetric properties of the ground state (and 
excitations) in different phases. In the context of our discussion a typical 
question from the first group looks like the following: what does tempera- 
ture deconfinement phase transition correspond to in terms of correlators? 
The second group provides questions like: where is screening of zero A^-ality 
charges at large distances hidden in the expression (jH))? We have no possi- 
bility to discuss these important issues in the present review and refer the 
reader to original literature and references therein (see review 

3 Mechanism of confinement and dual Meiss- 
ner effect 

3.1 Effective fields definition 

The formalism considered so far allows one to perform the expansion of Wil- 
son loop (jHl) and static potential (fT^ over the full set of field correlators (0), 
dZj), (fT^ . (fTTj) in the whole range of distances. In what follows we will use 
these results to calculate the effective confining field in hadrons and study 
some of its phenomenological applications'*. 

It is well-known that the static potential at small quark-antiquark dis- 
tances r -C ^QCD in Born approximation of perturbation theory has the 
form 



where = 4/3 is the quadratic Casimir operator in fundamental represen- 
tation. The color factor Ci;' is the only difference between this potential and 

"^Dynamics of effective fields is considered in [5^], [SZ| 




(18) 



r 
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Coulomb one in electrodynamics. One can introduce the field 



which has the meaning of the force acting on the quark. 
Let us define the effective field as follows, 

J^'^,{x) = {TTW{C))-\TTzgG,,{x,xo)W{C)). (20) 

Index J stresses that the field T'jj^^ix) is the functional of the external current 
J corresponding to W-loop W[C\ It will be demonstrated in the next section 
that this effective field at small distances is reduced to color-Coulomb field 

(Hi), (HI- 

Notice that one can write down the effective field using the connected 
probe Elj {TtW{C,Cp)), where 

WiC, Cp) = W{Cp, xMx, xoMxo, z)W{C, z)^z, xo)$(xo, x) (21) 

is the W-loop with the contour consisting of the (small) probe contour Cp 
connected with the contour C along some trajectory going through the point 
Xq. This quantity depends on the position of the "reference point" Xq as well 
as on the shape of the trajectory connecting C and Cp. We will choose the 
trajectory going along the shortest path from point x to the minimal surface 
of the W-loop, see Fig. 1. 

The effective field in the case of probe contour Cp with the infinitesimal 
surface 6a^i, can be written as 

J^i^ix) 5cr^,ix) = (It WiC))-' ((Tr WiC Cp)) - (Tr WiC))) ^ M(C, Cp) 

(22) 

In particular, if the probe contour has a size ax a, the relation for the electric 
field follows, 

„.5^(.).^(^^. (23) 

where n is the unit vector defining the orientation of probe contour in coor- 
dinate space^. 



^Let us remind in this context the expression for the moment of forces acting on the 
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3.2 Definition of effective currents 

In abelian theory Gpi^{x,xo) = F^y^x) and equation (pH) defines field dis- 
tribution satisfying Maxwell equations with the external electric current 

g''J,{x)=g^ j^dzp5^*){z-x) 

^^/,(^) = (24) 

where g denotes the electric charge. Let us now proceed with the nonabelian 
case. Using the differential relations for phase factors (see, e.g., [H] - |S21), 
one can formally write down the effective "electric" and "magnetic" currents 
as 

3-i{x) = (Tr W^(C))-i i^{i:i^x^;x)tgDpF^,{x)^{x-x^)W{C)) 

V rf,^!ff^^!f|^(Tr [G.,(«, Xo)G,.(a:, x,)]W{C))^ (25) 
ki{x)=g^{i:iW{C))-'^ 
jy^ du^Mdu,{s,x) ^^^ [G.,(«,Xo)G,.(x,a:o)] W{C)) (26) 

where the integration contour is given by the function m^(s,x) with the 
boundary conditions ^'^(O and square brackets denote 

commutators in color space. Index J indicates that the "electric" current 
and " magnetic" one fc^ are functionals of the external current J given by the 
W-loop. The currents defined in such a way can now be considered as sources 
of effective "electric" and "magnetic" fields according to effective "Maxwell 
equations" 

^ w^^i(^) = ; ^-^/m(^)=^.'(^)' (27) 

frame with the electric current / in the magnetic field B, known from the general physics. 
Namely, when the frame is oriented in the plane (n*^^), n^^^) and n^^^ is chosen orthogonal 
to magnetic field, the moment of acting forces M takes the form n^^^ ■ B — where 
3 = 1^. Comparing this relation with |SJ one sees that M{C, Cp) defined in (|22() means 
the "dual" moment of acting forces. 
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Equations and (j^HI) define effective currents wliicli satisfy (P7j) witfi 
tlie definition (j2(Jj) identically. Notice that in (j26|) nonabelian Bianclii iden- 
tities D^F^y{x) = respecting the gauge nature of QCD are used. It is 
obvious from ()27|) that both electric and magnetic effective currents are con- 
served since tensor JF^,^ is antisymmetric. 

Let us note that the lowest term of the W-loop expansion, which con- 
tributes to fc^(a;) (j26p . is proportional to nonabelian field strength correlator 
of third order. Therefore the value of magnetic current is proportional to 
correlator {EfB^E^)f°-^'^eijk, i.e. the effective magnetic current emerges due 
to the nonabelian emittence of the colormagnetic field by the colorelectric 
one The averages of the type (Tr Gai3{x,Xo)G^s{y,Xo)W{C)) in rhs of 
()25|) . ()26|) define the nonlocal gluon condensate in the presence of W-loop, 
which saturates to the constant value far from the W-loop. We do not ad- 
dress here an interesting question about possible microscopic nature of the 
currents (j^ . ^I^i . in particular, the question to what extent the magnetic 
current may be understood as corresponding to some propagating point- 
like particles, "abelian monopoles". Instead, we take (|^ . (pUj) as primary 
effective definitions. 

If the gauge coupling is small, one can use for the electric current 
the equation of classical gluodynamics. 



where J^{x) = J^(x)T", J^(x) = j^dz^5^^\z — x). In the leading order 
in gauge coupling = / (47r) the second term of does not contribute, 
and the expression for the electric current reads as 



i.e. it has a form of classical current of electrodynamics with the charge 
CpO-s- In particular case of static quark and antiquark Maxwell equation 
with this current reproduces the color-Coulomb potential (fT5|l . 




(28) 



(29) 
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3.3 Effective fields distribution in two-point approxi- 
mation 

Let us consider the rectangular W-loop of static quark and antiquark. Rely- 
ing on the hypothesis of bilocal (gaussian) dominance we take into account of 
only the bilocal correlator contribution to the effective fields assuming that 
higher correlators do not lead to essential modification of the confinement 
picture. The effective field in bilocal approximation reads 

T^,^{x)= / do^p{y)g^D^^l [x-y), (30) 
Js 

where y & S, S is the minimal surface of the W-loop, and bilocal correlator 
is defined in flT3|). 

Let us denote n = R/i? the unit vector directed from quark to antiquark 
and rewrite ()30p in the form 

:F^u{x)= [ d%Ti {gF^,ixMx,y)ngE{y)^y,x)) , (31) 
Js 

which clearly indicates that the magnetic field B is absent. The substitu- 
tion of parametrization (fT3j) for (jHTj) yields the following expression for the 
effective electric field, 

R oo 

Si{T, R) = n, jdl I dt (^6,,D{z) + , (32) 

-oo 

where z = {r — nl, t). The perturbative part of the field corresponding to the 
contribution of the formfactor Di to ((221) '^^^ be represented as the difference 



where 5^°"'(r) is the color-Coulomb field (fTH |) . (fT!7 |) . 



£^°"'(r) = (34) 



r 



The corresponding formfactor. 



Oriz) = (35) 
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can also be calculated directly in perturbation theory pUj . 

It was discussed in previous chapter that the confinement is the conse- 
quence of stochastic nature of gluon field fiuctuations, which reveal them- 
selves at separations of the order of the correlation length A and lead to the 
exponential fall off of the field correlators, see (|T7jl . One can show that since 
at large separations the string acts on quark with the force a, the formfactor 
D should be normalized according to 

«(^=)-^-p(-^)- (36) 

On substituting for one calculates the corresponding field, 

R/\ 













n — I dl 








) 


71" J 









(37) 



where Ki is the McDonald function. The string tension a can be considered 
as a scale QCD parameter (it is related to Aqcd through equation (j3^ ). 
Numerical value a ~ 0.18 GeV^ is determined phenomeno logically from the 
slope of the meson Regge trajectory, see e.g. [^. It is easy to verify that 
if the point x is placed at the symmetry axes, the relation between the field 
£^ and the nonperturbative part of the static potential corresponding to 
formfactor D Q, (fT^. which we denote , reads 

5^(0,R) = VV^(/?). (38) 



The distribution of the field \S{xi, 0,^3) | (jn2j) is shown in Fig. 2 at QQ- 
separation 2 fm. One can see at the figure the peaks of the color-Coulomb 
field ()34|1 over the quark and antiquark, and the string ()37|1 between them, 
with the universal profile £{p), 

Sip) = 2a{l + exp (-0 , (39) 

where p is the distance to the QQ axis. 
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3.4 Magnetic currents distribution and Londons equa- 
tion 



To perform more detailed analysis of the magnetic currents distribution 
in the case of static quark and antiquark let us apply the first Maxwell 
equation ^T7\ to the electric field in bilocal approximation (|?7j). 
Then one can seethat the magnetic current k has a form 



while the magnetic charge is absent. It is clear that the perturbative color- 
Coulomb field (fTUj) does not contribute to pUj) . A nonperturbative field (|H7jl 
is directed along the quark-antiquark axis, therefore magnetic current winds 
around the axis. In particular case of the saturated string the polar 
component of the magnetic current k^p takes a form 



One can see that the value of the current rises linearly near the axis and falls 
exponentially at large distances from it. 

The vector distribution of magnetic currents in the case of QQ-separation 
i? = 2 fm is shown in Fig. 3. This distribution resembles the one of the elec- 
tric superconducting currents around the Abrikosov string in superconduc- 
tors [91 , and is another hint in favor of dual superconductivity mechanism 
of confinement jBHl- An exponential behavior of current and field at large 
distances means that the dual Londons equation 



is satisfied. Indeed, the only component of the polar vector pij) is directed 
along z axes and has a form 



k = rot £, 



(40) 




(41) 



rot k = \'^£ 



(42) 



(rot k).(p) = i^=7(p)A-2£:(p), 




(43) 



where the universal profile S{p) is defined in (I39j) . and function 




(44) 
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rises monotonically from —2 and at p ^ A tends to unity as 7(p) ~ 1 — 3A/ p. 

One concludes that the confinement mechanism is related to cyclic mag- 
netic currents ()41|) squeezing the electric field into the tube of string with 
the exponential fall off outside it, and satisfying the dual Londons equation 

(m, (m. 

3.5 Vacuum polarization and screening of the coupling 
constant 

We turn now to the second Maxwell equation for the static quark and anti- 
quark, the Gauss law 

div S = p, (45) 
where the field £ ((221) is the sum 

and £:^i'°s<', £^ are defined in (jSHD, (jnil),(j27|), while the nonperturbative field 
£Di,np QYi?,uies the exponential fall-off of the formfactor Di at large distances. 
In this section we introduce additional assumption about charge distribution. 
This assumption is confirmed a posteriori by the lattice results in abelian pro- 
jected gauge theory (compare e.g. the distributions of effective field and its 
nonperturbative part along QQ axis in Figs. 6,7 with corresponding distri- 
butions in Fig. 21 from paper ^^). We assume that the nonperturbative 
contributions to the charge density cancel, 

divf:^i'°P = -div5^, ^47) 

so that the charge density has a form 

p = 47rCi.a,(5(r) -5(r-R)). (48) 

Using explicit expression for the field £^ dHIj), we find from (gTj) the ' screen- 
ing" charge density p(r), 

div£:^^'°P = p(r) -p(|r-R|), (49) 
*-) = -^rA-,(-). (50) 
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Relying on (gHl), dSHD, one calculates the field ''^p, 

5i?i,np ^ _ g(|r-R|)(r-R) 



r3 |r-R|3 ' 



where Q is the "screening" charge, 



Q{r) = =^ / x''Ki{x)dx, (52) 







It is obvious that there is no field at large distances from quark and antiquark 
due to confinement, and the full charge Q{r) defined as 

Q{r) = CFas{r)-Q{r), (53) 

turns to zero. The condition Q{r)\^_^^ = leads to the relation jiU] 

Cptts = SaA^ (54) 

between the strong coupling at large distances and parameters a, A respon- 
sible for confinement. The behavior of the charge Q{r) at standard values 
a = 0.18 GeV^, A = 0.2 fm and constant value = 0.42 calculated from 
fl54|) is shown in Fig. 4. The mean radius of the screening according to the 
figure is of the order of 0.5 fm. 

The behavior of the strong coupling taking into account the background 
confining fields was studied in [HHl EHI in the framework of the background 
perturbation theory. It was shown that the confining background leads to the 
modification of the logarithmic running of coupling according to as(g^) — >■ 
as{q'^ + where the "background mass" ms ^ I GeV~ A~^ is related to 
the energy of the valence gluon excitation and the gluon correlation length. 
One can see that at large distances r ^ A the background coupling tends 
to constant ("freezes"), while at small ones it turns to the running coupling 
of the ordinary perturbation theory. The relation between parameters now 
takes the form 

CpasiX) = 3(7 A^ (55) 

where as(X) is the freezing value of background coupling equal to the ordi- 
nary running coupling at scale A. This equation relates two alternative scale 
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parameters of quantum theory, Aqcd and the string tension a^. 

In Fig. 5 the background running couphng is shown by dotted curve. The 
behavior of the running charge Q'^'^^lr) is shown by sohd curve. As one can 
see from the figure, the effective charge has a maximum at r ^ 0.3 fm. 

Using standard values a = 0.18 GeV^, A = 0.2 fm and constant value 
as = 0.42 we plot the following field distributions. In Fig. 6 the projections 
of fields £^{0, 0, X3), 5(0, 0, 0:3) and £^^'"^%0, 0, X3) into the quark- ant iquark 
axis are shown. Note that the fields £^1'°p and £^^'°^'^ cancel in the middle 
of the string. In Fig. 7 the projections of fields £'^(0,0,X3), £^{0,0, X3) + 
£^^'"^(0, 0, X3) and £(0,0,X3) onto the quark-antiquark axis are plotted. In 
Fig. 8 the vector distribution of the displacement field £(xi, 0,2:3) is shown 
demonstrating that the field is squeezed in tube with the width of the order 
of A. In Fig. 9 the vector distribution of the solenoid field £'^(a;i, 0, 0:3) + 
f:^i'''P(xi,0,X3) is plotted. 

It is convenient to define the isotropic dielectric function e{r), 

e{r) = (56) 
One can calculate that at large distances r ^ A it is exponentially small. 



indicating the disappearence of the color-Coulomb field both inside and out- 
side the string. 



4 Hadrons with three static sources 
4.1 Green functions and W-loops 

Physical hadrons are nonlocal extended objects, therefore to construct their 
Green functions one should use nonlocal quark and gluon operators 

q\x,Y ) = q^{x)mx,Y), (58) 

^It was shown in that the value Aqcd = 241 MeV computed in lattice [Hll in MS 
regularization scheme with Uf — corresponds to the freezing value as (A) = 0.42, the 
latter being in complete agreement with H55|l . 
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gaix, Y) = gb{x)^abix, Y), (59) 

as well as the local one Gi{x) = ga{x)tf'^'' . Here and in what follows j, ... = 
1, 2, 3 are color indexes in fundamental representation, and a,b, ... = 1, .., 8 in 
adjoint one; ga denotes the valence gluon operator of background perturba- 
tion theory [^1111^], and Gl{x) transforms as Gl U^-'G-'-,Ul' under gauge 
transformations. One can construct gauge invariant combinations of these 
operators using symmetric tensors 6j, 6""^, d"'^'^ and antisymmetric ones Cijk, 

jabc 

BY{x,y,z,Y) = e,,kq\x,Y)q^{y,Y)q\z,Y), (60) 

G^4\x,y,z,Y) = r''g,{x,Y)g,{y,Y)g,{z,Y), (61) 

G^^\x,y,z,Y) = d'''^ga{x,Y)g,iy,Y)g,{z,Y), (62) 

GA{x,y,z) = Gi{x)^';{x,y)Gi{y)^riy,^)G';.iz)Kiz,x). (63) 

First three constructions have a structure of F-type with the string junction 
at point Y, where the color indexes are contracted with the (anti-) symmetric 
tensor, and the latter one has a structure of triangular type. Let us stress 
that the wave function of triangular type is possible only for glueballs but 
not for baryons, see jHIj. 

Hadron Green function has a form 

^,(X,X) = (vI/+(X)vI/,(X)), (64) 

where = G^, Gy, By] X = x,y, z in the case of Ga and x, y, z, Y for Y- 
states. The vacuum average (...) leads to the product of Green functions of 
quarks or valence gluons, which are proportional to the parallel transporters, 

{qj{x)q\x)) ~ $}(x,x), 

{ga{x)gb{x)) ~ ^ab{x,x). (65) 

Therefore the hadron Green function is proportional to the W-loop of this 
hadron, see equation (fTT|) . and is reduced to the W-loop in the case of static 
sources. W-loops of baryon and F-type glueball take forms correspondingly 

= ^(e,,fce^'^''^'«fMC^i)$MC'2)'^'''(C^3)), (66) 



24 



Wa''^ = ^(rt"'''''='$"°'(Ci)$*^'(C2)$'^'='(C3)). (68) 

Trajectories Ci formed by the sources are shown in Fig. 10. A W-loop of 
A-type glueball at large distances can be represented as a product of three 
meson W-loops 97], 

W^(X, X) = W{x, y\x, y)W{y, z\y, z)W{z, x\z, x). (69) 

Corresponding contours are shown in Fig. 11. 



4.2 Static potentials 

Static potentials of hadrons with three static sources are calculated in bilocal 
approximation of the field correlator method [HZl I103j in the same way as 
meson ones^. For hadrons of F-type let us denote n^"^ the unit vector directed 
from the string junction to the a-th quark and Ra the separation between 
this quark and the string junction. The potential in baryon reads 

(\ pRa rRb /"OO 

V-V nf^nf / / dldl' dm,,,{z,,), 
a=b a<bj Jo Jo Jo 

(70) 

where Zab = {In'-"^ - l'n^^\t). One can represent this potential in the form 

Vs = V' + V'^ + V"^, (71) 
where is the color-Coulomb potential 

Tij is the 2-th and j-th quark separation. We will take into account the charge 
screening replacing in (f?^ Cpois with Q defined in fISHj) . (jH21)- Terms and 
V^'^'^ denote the diagonal and nondiagonal parts of the potential corresponding 
to the correlator D; \s determined by the first and I^"*^ by the second sum 

^The effect of cliarge screening is taken into account in (HHI 
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in (j7(Jj) . One can find explicit expressions for and V^'^ in |l(J3j . We just 
note liere tliat is a sum of quark-antiquark potentials ©, (HSl), 

y'^(/?l,i?2,i?3) = $^^''(^a). (73) 

a 

Characteristic feature of the potential (fTOj) is an increase of its slope when 
the source separations are increasing. In Fig. 12 the behavior of the baryon 
potential with the color-Coulomb part subtracted is shown in comparison 
with the lattice data |lU4j as a function of the total length of baryon string 
L = Ra- A tangent with the slope a is shown by points. One can see 
from the figure that the potential slope becomes significantly less than a at 
L < 1 fm. This effect is induced by the influence of the correlation length 
of confining fields |lU3j . In Fig. 13 a dependence of the baryon potential in 
equilateral triangle on the quark separation is given in comparison with the 
lattice data |lU5j . Note the agreement between analytic and lattice calcula- 
tions within the accuracy of a few tens MeV. For the F-glueball potential 
the Casimir scaling holds, 

YL = 9l (74) 

where C3 = {N^ — l)/2Nc = Cp and Cs = Nc are quadratic Casimir operators 
in fundamental and adjoint representations. 

A potential in A-glueball in the case of equilateral triangle with the side 
r has a form [HZj 

Viir) = ^V%r) + V'ir) - 2y°<^(r). (75) 

Let us note that and V^'^ depend on the valence gluons separation but 
not on the separation between the gluon and the center of the triangle, and 
that the term —2V^'^ corresponds to the interaction of three effective quark- 
antiquark W-loops. The behavior of the potentials Vq and in equilateral 
triangle in dependence on the source separation r is shown in Fig. 14. The 
potential Vq goes above the Vq because of the positive contribution of non- 
diagonal term V^'^ to Y-type glueball and negative to triangular one, as well 
as for the greater slope of diagonal term in the case of F-glueball. 
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4.3 Fields distributions 



The field in baryon is defined jlO] as the square average 

{S^''^? = \{{S^)f + {Sl^f + (76) 
of fields calculated for the probe plaquette joint to the trajectory Cj, 

£fi)(x, R(i), R(2), R(3)) = £:^^(x, R(i))-^ £^-^(x, R(2))-l £^^(x, R^^)). (77) 

Normalizing coefficient 2/3 in ()76|) is chosen due to the condition that at large 
separations the field acting on quarks equals to a. According to (f76|) . (ffTjl . 
the field in baryon is expressed through fields of effective quark- ant iquark 
pairs, with positions of antiquarks coinciding with the string junction. The 
distribution of the field E^^"^ taking into account only the contribution of 
formfactor D is shown in Figs. 15 and 16 in the plane of quarks forming an 
equilateral triangle with the side 1 fm and 3.5 fm respectively. One can see in 
Fig. 16 three plateau with the saturated profile, and small growth of the field 
around the string junction point, the relative difference of values amounts to 
1/16. A surface formed by the confining field with the value a is shown in 
Fig. 17 for quark separations 1 fm. One can see the small convexity in the 
region of the string junction. 

A field in A-type glueball is a sum of meson fields with gluon pairs acting 
as the effective sources jlUj . 

Sf\^, r(i), r(2), r(3)) = ^ S'\^ - r^, r^^+^^^^d^ _ ^W)^ (73) 

i=l 

where r^') denotes the position of i-th valence gluon. In Fig. 18 the field 
distribution |£^''(x)| in valence gluon plane is shown at gluon separations 1 
fm, and in Fig. 19 the surface l^^i (^)l = cr is plotted for the same gluon 
separations. Let us note that according to one can calculate the static 
quark-antiquark potential as the work of the force acting on the quark, done 
on separating the latter from the antiquark to the distance R. Analogous 
relation are valid for the field and potential of A-type glueball, see (|75|l . (f78|l . 
It is clear that the nondiagonal part of the A-glueball potential V^'^ equals to 
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the work of force acting on the effective quark from the external string and 
is therefore related to the interference of the meson fields in the vicinity 
of valence gluons of the order of A. 

5 Conclusions 

In this paper we have systematically treated: the vacuum fields in QCD, 
the confinement mechanism, the QCD string formation and finally, the field 
distribution inside hadrons. Everywhere we have used the field correlators as 
a universal gauge-invariant formalism, which allows to describe all phenom- 
ena appearing in QCD. In description of vacuum fields the most important 
property is the Gaussian dominance: the lowest (Gaussian) correlator is dom- 
inating on the minimal area surface of the W-loop, and there are sufficient 
grounds for the statement that the total distribution of higher correlators 
does not exceeed few percent. This phenomenon, found on the lattice 
is not yet fully understood, see [TTl HH], although it gives an explicit dy- 
namical picture, which possibly is incompatible with the old physics of the 
instanton gas, of Z2-fiuxes etc.. Therefore, one can assert that the picture of 
the maximally stochastic QCD vacuum is a very good approximation to the 
reality. One can remember that the measure of coherence is associated with 
the weight of the contribution of higher correlators, e.g. for the instanton 
gas the total contribution of higher (non- Gaussian) correlators is dominating. 
Moreover, the vacuum correlation length A (i.e. the factor in the exponent 
for the asymptotics of the Gaussian correlator) is relatively small, A ~ 0.2 
fm for the quenched vacuum. This value is much smaller than the typical 
hadron radius, ~ 1 fm. Theoretically, the smallness of A is connected to a 
large mass gap for glueballs and gluelumps, since A = 1/M, where M is the 
lowest gluelump mass, M ~ 1.4 GeV, which is calculated both analytically 
and on the lattice [T^OITi] . 

Let us turn now to the confinement mechanism. From the point of view of 
field correlators, confinement occurs due to the appearence of a specific term 
in the Gaussian correlator, denoted D(a;^), which violates Bianchi identities 
in the abelian case and therefore is absent in case of QED. If, however, one 
considers the U{1) theory with magnetic monopoles present in the vacuum. 
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then the function D{x'^) is nonzero and is proportional to the monopole cur- 
rent correlator. The next step is to find the source of D{x'^) (i.e. the source of 
confinement) in the nonabelian theory. It was done in |lU6t I1U7] . where the 
derivatives of D(x^) were connected to the triple correlator (EEB). Thus 
the problem of establishing of the confinement mechanism in the formalism 
of field correlators reduces to the problem of calculating D{x^) and the triple 
correlator and to the finding the conditions of its appearence/disappearence 
(e.g. in QCD - as functions of temperature or baryon density). The lat- 
tice calculations confirm the disappearence of -D(x^) at the deconfinement 
temperature T^, and with it have confirmed all cofinement picture in the 
framework of the present method. One expects that at the next step- by com- 
puting correlators (including D{x'^)) with the help of the gluelump Green's 
functions in the whole x region - one will make the field correlator method 
self consistent, and the problem of confinement will be solved quantitatively 
and in principle. 

At the same time, this universal formalism of field correlators can be used 
to study the distribution of effective fields and currents, defined with the help 
of the W-loop. This representation, see chapter 3, enables one to describe, on 
one hand, the dual Meissner effect [1111111, and on the other hand, it relates 
to the effective Lagrangian approach of Adler and Piran [22] and dielectric 
vacuum models, P7| |2H] and subsequent papers. Indeed, the field correlator 
method not only admits this approximate qualitative interpretation, but also 
yields explicit expressions for the density of effective electric charges and ef- 
fective magnetic currents. Being the gradient of the color-Coulomb potential 
at small distances, the effective field condenses into a tube on the characteris- 
tic hadron scale and ensures confinement. In the process the strong coupling 
constant is screened due to the vacuum polarization by nonabelian gluonic 
interactions. 

Finally, let us summarize the contents of the last chapter devoted to 
field distributions inside hadrons with three constituents (sources). Here 
the field correlator method is the only quantitative analytic method, and its 
comparison with numerical (lattice) results is very interesting. One can note 
that in the method one has only two parameters - the string tension a and 
the correlation length A, A being expressed through a, and a is playing the 
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role of the scale parameter related to Aqcd- The baryon potential computed 
in this way jl(J3j is in good agreement with lattice calculations and gives an 
independent confirmation that baryon strings have the structure of the Y- 
type with the string junction. Moreover, the field correlator method explains 
the smaller slope of the baryon potential at the typical hadron distances, 
known from the baryon phenomenology - the decrease of the slope is caused 
by the string interference effects connected to nonzero correlation length A. 
The three-gluon glueballs, in contrast to baryons, can have the structure of 
both y-type and A-type ^7\. However, the latter is preferred energetically. 
In the concluding part of the last chapter the field distributions in baryons 
and in the A-type glueballs are given, where one can visualize the shape of 
the string in these hadrons. 

Summarizing, one can say that the universal language of the field corre- 
lator method turns out to be extremely convenient in all cases considered. 
In particular, it enables one to formulate the gauge-invariant description of 
the QCD vacuum as some medium with properties ensuring confinement. 
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Q 

Figure 1: A connected probe (PT|) for static quark and antiquark 




Figure 2: A distribution of the field |5(xi,0,X3)| (j32|l at quark- antiquark 
separation 2 fm. Gutted peaks of color- Coulomb field and string between 
quark and antiquark are clearly distinguished. The standard values of pa- 
rameters a = 0.18 GeV^, A = 0.2 fm are used. 
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Figure 3: A vector distribution of magnetic currents (|H7j). pUj) at quark- 
antiquark separation 2 fm. Positions of quark and antiquark are shown by 
points. 
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Figure 5: A running background coupling CpCiBir) (dotted curve) and 
running effective charge Q'^'^^ = CpOiBij) — Q{r-) (solid curve) vs. the dis- 
tancefrom the quark. 
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Figure 6: Distributions of the projections of the fields £^{0,0,X3) (sohd 
curve), £^i'°P(0,0,a;3) (dashed curve) and £^i'°s^(0, 0, X3) (dotted curve) 
onto the quark-antiquark axis at QQ-separation 3 fm. 
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Figure 7: Distributions of the projections of the fields £ (0,0,2:3) (solid 
curve , £^{0,0,X3) + £^i'"P(0,0,a;3) (dashed curve) and £{0,0, X3) (dotted 
curve) into the quark-antiquark axis at QQ-separation 3 fm. 
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Figure 8: Vector distribution of the field S{x\, 0, X3). Positions of quark and 
antiquark are marked by points. 
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Figure 9: Vector distribution of the solenoid field S^{x\,^,xz) + 
0, X3). Positions of quark and antiquark are marked by points. 
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Figure 11: A W-loop of A-type. 
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Figure 12: A potential in baryon (f7n|) with the color-Coulomb part contracted 
(solid curve) in comparison with the lattice data |lU4j points in dependence 
on the total length of the baryon string L. A value of the string tension is 
a = 0.22 GeV^. According to ()55|) . the corresponding value of correlation 
length is A = 0.18 fm. 
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Figure 13: A dependence of the baryon potential in equilateral triangle on 
quark separation r (solid curve) in comparison with the lattice data |lU5j 
(points). A value of the string tension is a = 0.17 GeV^ (A = 0.21 fm). 
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(dotte curve) in equilateral triangle vs. the sources separation r. 
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Figure 15: A distribution of the field £^^^ dZEI), dZZI) in GeV/fm with the 
only correlator D contribution considered in the quark plane for equilateral 
triangle with the side 1 fm. Coordinates are given in fm, positions of quarks 
are marked by points. 
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Figure 17: A surface (x)| = cr at quark separations 1 fm. Coordinates 
are given in fm, positions of quarks are marked by points. 
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Figure 18: A distribution of the field \S\^ (x) | (f75)) in GeV/ fm of the triangu- 
lar glueball in the plane of valence gluons with separations 1 fm. Coordinates 
are given in fm, positions of valence gluons are marked by points. 
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Figure 19: A surface (x)| = cr at valence gluons separations 1 fm. 
Coordinates are given in fm, positions of valence gluons are marked by points. 
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